Spin Hall transistor with electrical spin injection 
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The realization of a viable semiconductor transistor and information process- 
ing devices based on the electron spin has fueled intense basic research of three 
key elements: injection, detection, and manipulation of spins in the semicon- 
ductor channel. The inverse spin Hall effect (iSHE) detection of spins injected 
optically in a 2D GaAs^ and manipulated by a gate- voltage dependent internal 
spin-orbit field has recently led to the experimental demonstration of a spin 
transistor logic device.*^ The aim of the work presented here is to demonstrate 
in one device the iSHE detection combined with an electrical spin injection and 
manipulation. We use a 3D GaAs channel for which efficient electrical spin injec- 
tion from Fe Schottky contacts has been demonstrated in previous works.l^^ In 
order to experimentally separate the strong ordinary Hall effect signal from the 
iSHE in the semiconductor channel we developed epitaxial ultrathin-Fe/GaAs 
contacts allowing for Hanle spin-precession measurements in applied in-plane 
magnetic fields. Electrical injection and detection is combined in our transistor 
structure with electrically manipulated spin distribution and spin current which, 
unlike the previously utilized electrical manipulations of the spin-orbit field^ or 
ballistic spin transit time,^^ is well suited for the diffusive 3D GaAs spin chan- 
nel. The magnitudes and external field dependencies of the measured signals 
are quantitatively analyzed using simultaneous spin detection by the non-local 
spin valve effect!^^ and modeled by solving the drift-diffusion!^^ and Hall-cross 
response equations for the parameters of the studied microstructure. 

The pioneering works on electrical spin injection and detectiori^^*^ in non-magnetic chan- 
nels were done in metals, taking advantage of the compatibility with conventional metal 
ferromagnets which formed the injection and detection electrodes. The non-local spin valve 
effect, utilized in these seminal studies, measures the dependence of the electro-chemical 
potential at the detection ferromagnetic electrode on the relative orientation of the magne- 
tization of the electrode and the spins in the non-magnetic metal underneath it. 

A fundamentally distinct method for the electrical detection of spin currents in non- 
magnetic conducting channels is based on the iSHE. The approach was first demonstrated in 
metal devices with out-of-plane magnetized ferromagnetic injection contacts and compared 
in the same microstructure with the non-local spin valve signal .'^^'^ The iSHE detection does 
not utilize a reference ferromagnetic probe. Instead, a transverse spin dependent voltage 
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of a Hall cross fabricated directly in the non-magnetic channel provides the measure of the 
out-of-plane spin polarization of the propagating electrons. The basic physics distinction 
between the two approaches is that the spin valve effect originates from the exchange- 
splitting of carrier bands in the ferromagnetic probe while the iSHE originates from the 
spin-orbit coupling in the non-magnetic conduct or .f^Slf^ni 

Several recent works have applied the non-local spin valve method for detecting spins 
in the GaAs semiconductor."^ We have prepared specially designed Fe/GaAs microstruc- 
tures which allow us to reproduce these previous non-local spin valve experiments and to 
simultaneously demonstrate in the same microdevice the detection of the spin current by 
the iSHE in the lateral GaAs microchannel. An important feature of our Fe/GaAs devices 
which enables the iSHE detection is the strong in-plane anisotropy of the ultrathin-film 
(2 nm) Fe electrodes. It allows us to apply sufficiently strong in-plane magnetic fields along 
the Fe hard- axis for performing Hanle spin precession experiments, without aligning the Fe 
magnetization with the external field. This is the suitable geometry for the iSHE detection 
because spins injected from the in-plane magnetized Fe electrode precess in GaAs out of 
the plane of the transport channel. Furthermore, the iSHE and the ordinary (Lorentz force) 
Hall effect contributions can be experimentally separated in this set up which is essential for 
detecting the iSHE in semiconductors. Recall that only in the high carrier density metals, 
the ordinary Hall effect is relatively weak and can be neglected in the iSHE experiments in 
external magnetic fields.^ 

Another important feature of our device is that it allows to manipulate the distribution 
of the diffusive spin current by applying electrical bias across the transport channel. Our 
work complements the previous realization of the GaAs spin transistoi^^ in several aspects, 
(i) While the device reported in Ref. 2 acted as a spin analogue of a field-effect-transistor 
logic element, the present device is rather a spin amplifier, (ii) Our device does not utilize 
optical spin injection but electrical injection from a ferromagnetic contact, (iii) We measure 
simultaneously the output spin current by the iSHE and the output spin polarization by 
the non-local spin valve effect which allows us to quantitatively analyze the detected spin 
signals. 

Our Fe/n-GaAs heterostructure was grown epitaxially in a single molecular-beam-epitaxy 
chamber without breaking ultra high vacuum conditions during the whole growth process. 
The sample was deposited on an undoped GaAs substrate. The heterostructure comprises 
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250 nm of low Si-doped GaAs (5x 10^^ cm~^), followed by 15 nm of GaAs with graded doping, 
and 15 nm of highly Si-doped GaAs (5 x 10^^ cm~^). The doping profile yields a narrow 
tunnel Schottky barrier between GaAs and Fe favourable for spin injection/detection.'^'^ 
The growth temperature of GaAs was 580° C. The sample was then cooled to 0°C for the 
growth of the 2 nm Fe layer. The refiection high energy electron diflFraction pattern observed 
after the Fe deposition confirmed the epitaxial growth of cubic Fe. The Fe film was capped 
by a 2 nm Al layer to prevent Fe oxidation. Electron-beam lithography and wet chemical and 
reactive ion etching were used to pattern the lateral GaAs channel with the Hall crosses and 
magnetic electrodes, shown schematically in Fig. la (for more details see Supplementary 
information). The distance between Fe electrodes is 4 /xm and between the Fe electrode 
and the Hall cross the distance is 2 /xm. Non-magnetic Au-contacts are located 100 /xm 
away from the central Fe electrode on each side of the channel. All measurements shown 
in this paper were performed at 4.2 K. The reproducibility of our experimental data was 
confirmed by performing measurements in three diflFerent samples with the same nominal 
heterostructure parameters and microdevice geometry and in each sample by contacting 
different combinations of available electrodes. 

The magnetic anisotropy of our Fe electrodes has the strong out-of-plane component 
(2 T) due to the thin-film shape anisotropy, the cubic magnetocrystalline component, and 
an additional uniaxial magnetocrystalline anisotropy originating from the broken [110] /[1 10] 
symmetry of our ultrathin Fe on GaAs. These anisotropics make the [110] in-plane crystal 
direction (y-axis in Fig. la) the easy magnetic axis with an anisotropy field of 0.2 T re- 
quired to align the magnetization with the [110] in-plane hard-axis (x-axis in Fig. la). This 
anisotropy field is significantly larger than external magnetic fields applied in the Hanle 
precession experiments with typical amplitudes up to 50 mT. When applied along the in- 
plane hard- axis, the component of the magnetization along the easy- axis is reduced by less 
than 10%. (For a more detailed discussion of the magnetic anisotropics see Supplementary 
information) . 

For characterizing individual Fe/GaAs Schottky contacts we performed three-point mea- 
surements between an individual Fe electrode and the two Au electrodes. These measure- 
ments confirm the tunneling nature of transport through the Schottky barrier. By measuring 
the tunneling anisotropic magnetoresistance of individual Fe/GaAs Schottky contacts, we 
inferred the strength of the above anisotropy fields and the technique also allows us to infer 
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Figure 1: (a) Schematics of the device for electrical spin-injection via a bias current and the 
simultaneous detection of the spin current by the iSHE voltage Vh and the spin accumulation by 
the non-local spin valve voltage Vnl- Green symbols show the notation corresponding to positive 
voltages and current, (b) Vnl measured by sweeping the external magnetic field By along the 
Fe easy-axis. The left (right) black arrow in each pair of arrows inside the plot indicates the 
magnetization state of the detection (injection) Fe electrode. The black and red arrows next to 
the label of the S^-axis indicate the field sweep direction. The data were measured at Ib = 
100 /J.A. (c) The difference between the non-local spin-valve signals at antiparallel and parallel 
configurations of magnetization in the injection and detection Fe electrode as a function of the 
spin-injection bias current, (d) Hanle spin precession/dephasing measurement by the non-local 
spin valve voltage in an applied in-plane hard-axis field Bx- As in panel (b), the arrows inside the 
plot indicate the magnetization states of the detection and injection Fe electrodes preset before the 
Hanle measurements. The color coding of the arrow pairs corresponds to the color of the measured 
data. The data were measured at Ib = 300 /iA and magnetic field sweep rate 2 mT/min. (e) 
Measured iSHE signal under same experimental conditions as in (d). The arrows notation is the 
same as in panel (d). 

the order in which the Fe electrodes switch their magnetization (see Supplementary infor- 



5 



mation). The magnetic characterization of each individual Fe electrode is essential for the 
interpretation of the spin injection, manipulation, and detection measurements which we 
now discuss in detail. 

In Fig. lb we plot the non-local spin-valve signal measured while sweeping the exter- 
nal magnetic field By along the easy axis of the Fe electrodes. In this geometry, the field 
triggers magnetization reversal in the Fe electrodes via a domain nucleation and propaga- 
tion process. The reversal fields are diflFerent in different Fe electrodes which allows us to 
control independently the magnetization orientations in the injection and detection elec- 
trodes. In the experiment, a bias current Ib driven between the injection (right) Fe/GaAs 
Schottky contact and the right Au electrode generates spin-accumulation underneath the 
spin-injection contact. A resulting diffusive spin current propagates into the unbiased part 
of the semiconductor channel with the Hall cross and the detection (central) Fe electrode. 
The lower value of the non-local voltage, measured between the detection Fe electrode and 
the left Au electrode, corresponds to parallel orientations of magnetization of the injection 
and detection Fe electrodes; the higher value of the non-local voltage corresponds to an- 
tiparallel magnetizations. These data together with the dependence of the amplitude of the 
non-local voltage on the bias current plotted in Fig. Ic, reproduce previous results of 
spin injection experiments in GaAs channels with Fe Schottky contactsJ ^ ^ ^ ^ ^ * ^ * Note that at 
Ib > 0, which is more favorable for efficient spin injection, electrons in the biased part of 
the channel drift in the direction from GaAs to Fe, i.e., are spin-selectively extracted from 
the semiconductor. 

Non-local spin valve measurements in magnetic fields Bx applied along the Fe in-plane 
hard axis are shown in Fig. Id. The black curve shown in the figure was obtained by setting 
the magnetizations in the Fe injection and detection electrodes in the parallel configuration 
before sweeping the in-plane hard- axis field; the blue curve was measured in the antiparallel 
magnetization configuration. At zero field we obtain the higher value of the non-local voltage 
for antiparallel magnetizations, consistent with Fig. lb. 

The Hanle dependence of the non-local spin valve signal on the hard-axis field B^ shown 
in Fig. Id refiects the precession and dephasing of spins in the GaAs channel, as quantified 
in detail in the theory section. The injected spins precess in the plane perpendicular to 
the applied field Bx^ i.e., acquire an out-of-plane component. Our observation of the iSHE 
signal due to the out-of-plane polarized spin current in the GaAs is demonstrated in Fig. le. 



6 



Consistent with the iSHE interpretation, the signal in Fig. le is zero at zero apphed field 
since in this case the in-plane polarized injected spins do not precess in the GaAs channel, 
i.e., do not acquire the out-of-plane component. The variations of the iSHE signal in Fig. le 
and of the Hanle non-local spin valve signal in Fig. Id occur at a comparable magnetic field 
scale. This confirms the precession origin of the out-of-plane spin component detected by 
the iSHE voltage. A full quantitative modeling of these Hanle non-local spin valve and iSHE 
curves is discussed in the theory section. 

The iSHE curves shown in Fig. le were obtained by subtracting the measured signals 
for oppositely preset polarizations of the Fe injection electrode before performing the Hanle 
measurement. The individual raw data measured for a given orientation of the injector po- 
larization have a linear contribution from the ordinary Hall eflFect. Importantly, the ordinary 
Hall contribution is independent of the preset orientation of the magnetization of the injec- 
tion electrode and is therefore experimentally removed by subtracting the data for opposite 
injector polarizations. The signals shown in Fig. le are, therefore, of pure spin origin and 
are due to the iSHE which is an odd function of the polarization of injected electrons. (For 
the comparison of the measured raw data and the iSHE signal and the discussion of the 
ordinary Hall eflFect in our sample geometry see Supplementary information). Since spin 
detection by the iSHE is performed directly in the GaAs channel, the corresponding signal 
depends on the magnetization state of the Fe injection electrode and, unlike the non-local 
spin valve signals, it is independent of the state of the Fe detection electrode. This expected 
behavior of the iSHE data is confirmed in Fig. le which shows measurements for diflFerent 
magnetization configurations of the injection and detection electrodes. 

For clarity, the Hanle non-local spin valve data shown in Fig. Id are symmetrized and, 
similarly, the iSHE data in Fig. le are antisymmetrized with respected to = 0. In 
Figs. 2b, c we show that the symmetrization (antisymmetrization) does not affect the key 
features of the measurements (see also Supplementary information). Nevertheless, the abso- 
lute value of the raw data have a noticable asymmetry which we attribute to the presence of 
dynamic nuclear spin polarization effects. 1^ To highlight the presence of the nuclear spins 
we plot in Figs. 2d,e the Hanle curves measured at different magnetic field sweep rates. For 
the applied in-plane field the Hanle curves show strong dependence on the sweep rate, 
consistent with the presence of the nuclear spin effects. On the other hand, the Hanle curves 
measured in the out-of-plane field B^ are independent of the sweep rate, implying that the 
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Figure 2: (a) Schematics of the experimental setup, (b) As measured (black) and symmetrized 
(red) non-local spin valve signal in the Hanle experiment with the in-plane hard-axis field Bx- 
The black arrow next to the label of the S^^-axis indicates the field sweep direction. The data 
were measured at Ib = 300 and magnetic field sweep rate 3 mT/min. (c) As measured 
(black) and antisymmetrized (red) iSHE signals for the same experimental conditions as in (b). 
(d) Symmetrized non-local spin valve signal in the Hanle experiment with the out-of-plane hard-axis 
field Bz measured at different sweep rates. The solid line shows the solution of the drift-diffusion 
equations discussed in the theory section. (e),(f) Symmetrized non-local spin valve measurements 
and antisymmetrized iSHE measurements in the in-plane hard-axis field Bx at different field sweep 
rates. 

nuclear spins do not significantly contribute in this experimental geometry. 

The role of nuclear spins in the measured Hanle curves can be understood from the 
following expression for the Overhauser field they produce and which renormalizes the 
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total effective field acting on electron spins 

B„ = fbr, (b • (S)) B . (1) 

Here / < 1 is a nuclear-spin relaxation leakage factor, 6^ = —8.5 T in GaAs, B is the unit 
vector of the external magnetic field, and (S) is the mean electron spin polarization. Because 
of the large out-of-plane magnetic anisotropy of the thin-film Fe (with the corresponding 
2 T anisotropy field) , the applied field in the out-of-plane Hanle experiment is not strong 
enough to significantly tilt the magnetization in the injection electrode from the easy-axis (^- 
direction). The projection of the injected electron spins to the applied field remains small 
in this experimental geometry which explains the absence of nuclear spin effects in Fig. 2d. 
For the in-plane Hanle experiment, the projection is given approximately by {S)Bx/Ba^ 
where Ba = 0.2 T is the in-plane anisotropy field. By fitting the solution of the drift- 
diffusion equation for electron spins to the measured non-local spin valve signals, shown 
in Fig. 2d and discussed in detail in the theory section, we infer that {S) ^ 5 %. This 
combined with the experimentally determined sign of the injected electron spin polarization 
yields B^ ^ Bx^ i.e., the total effective field experienced by the electron spins in GaAs is, 
B^ff ^ 2Bx- (The sign of the injected electron spin polarization, which is opposite to the 
magnetization in the Fe injection electrode, was obtained using the tilted magnetic field 
Hanle measurements'" which are discussed in the Supplementary information.) Consistently, 
the field scale on which the slow sweep rate Hanle curve varies in Fig. 2e is about a factor 
of 2 smaller than the field scale in Fig. 2d. 

Figs. 2e,f show the correspondence between the non-local spin valve and iSHE signals 
measured at different sweep rates of Bx- The correspondence demonstrates that the Over- 
hauser field acts in both measurements which provides further evidence that the measured 
Hall signals are of spin origin. We conclude the discussion of nuclear spins by pointing out 
that at sufficiently low sweep rates the role of the nuclear spins is merely in rescaling the 
effective magnetic field acting on electron spins in GaAs. Apart from the rescaling, the 
nuclear spins do not obscure the results of our spin transport experiments. This applies to 
both the non-local spin valve and the iSHE detection measurements. 

We now inspect the key symmetry of the iSHE signal which is the change of the sign of the 
Hall voltage upon reversing the spin current in the GaAs channel. This sign change is already 
seen by comparing Figs, le and 2c. In Fig. 1 we show data obtained from measurements 
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in which the right Fe electrode is used for spin injection and the left Fe electrode for spin 
detection (see Fig. la). Measurements shown in Fig. 2 were performed with reversed roles 
of the two Fe electrodes and, therefore, with the reversed orientation of the spin current 
(see Fig. 2a). As expected, this has no effect on the sign of the measured non-local spin 
valve voltage while the iSHE voltage changes sign for the two spin-current orientations. In 
Fig. 3 we provide an additional consistency check of the sign of the iSHE voltage. Here we 
show measurements in which the sign of the spin current is reversed by using the same Fe 
electrode for injection but biasing it with the left or right Au contact, respectively. The 
corresponding experimental setups are shown in Figs. 3a, b and the data plotted in Fig. 3c 
confirm that the sign of the measured Hall voltages is opposite for opposite orientations of 
the spin current in the GaAs channel. 
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Figure 3: (a),(b) Schematics of the experimental setups with opposite orientations of the spin 
current in the GaAs channel, (c) iSHE measurements in the in-plane hard-axis field Bx showing 
opposite signs for the opposite spin-current orientations. The data were measured at Ib = 300 fiA 
and magnetic field sweep rate 3 mT/min. 

In experiments shown in Figs. 1-3, the spins in the GaAs are manipulated by the external 
magnetic field via the Hanle spin precession. The spin current in these measurements is 
purely diflFusive in the part of the GaAs channel between the injection and detection Fe 
electrodes. On the other side of the channel with the bias current both the diffusion 
and drift are present. In Fig. 4 we show measurements in which we apply a bias between 
the two Au electrodes causing an additional drift current component Id on both sides of the 
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injection electrode. The corresponding experimental setup is shown in Fig. 4a. The current 
Ib driven through the injection Fe electrode is kept at a constant value of 300 /a A while the 
additional current is set to and ±100 /xA. The experiments illustrate the possibility 
to manipulate spins in our GaAs channel electrically via the bias dependent drift. As seen 
in Figs. 4b, c, both the spin polarization measured underneath the Fe detection electrode 
and the spin current measured by the iSHE depend on the applied bias between the Au 
electrodes. 

Qualitatively, the experimental data shown in Figs. 4b, c can be explained by a shift of 
the injected spin polarization profile from the injection electrode in the direction towards the 
Fe detection electrode in the case of Id = +100 fiA. In the experiment with Ijj = —100 /xA, 
the drift acts against diffusion on both sides of the injection electrode which makes the 
spin polarization profile decay more rapidly as we move away from the injection point. 
This explains why the detected spin signals are enhanced for positive and suppressed 
for negative Id- Our experiments demonstrate a method for modulating the output spin 
signal by electrical means. Among conventional transistors we can therefore find a loose 
analogy with the bipolar transistor amplifier. The spin current (spin polarization) detected 
by the iSHE (non-local spin valve effect) is the spin counterpart of the collector current and 
the additional drift current Id is reminiscent of the base current in the bipolar transistor. 
There is, however, a significant basic physics difference between the charge and spin based 
device. The latter uses the property that spin is not conserved. By applying the drift to 
electrons, the non-uniform spin-polarization profile along the channel can be shifted and the 
corresponding spin current increased or decreased which causes the electrically controlled 
modulation of the output signal. Note that our electrical spin manipulation is physically 
distinct from the electrical bias effect utilized previously in the ballistic Si spin channel.'^ 
The experiment in Si relied on the long spin lifetime and was based on electrically controlling 
the electron transit time through the channel relative to the Hanle precession time in an 
external magnetic field. 

We now proceed to the quantitative theoretical discussion of our spin injection, manipu- 
lation, and detection experiments. The spin dynamics in the GaAs channel can be modeled 
by the spin drift-diffusion equations. For the applied in-plane hard- axis field the spins 
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precess in the y — z plane and the corresponding Hanle curves are obtained by solving, 
dsy{x) d / j^dsy{x) \ Sy{x) 



dt 



where the nuclear Overhauser field is included in the total eflFective field B^^^ ^ as described in 
the discussion of Fig. 2. Analogous equations apply for the Hanle curves in the out-of-plane 
field Bz. In Eqs. ([5]), is the diffusion constant, Vd is the drift velocity, is the spin- 
dephasing time, g is the Lande-factor of electrons in GaAs, and /i^ is the Bohr magneton. 
The right-hand side of Eq. ^ for the Sy component describes the rate of spins parallel to the 
Fe magnetic easy- axis (^-axis) injected from the Fe contact to the GaAs channel at x = 0. 

In our experiments, the drift velocity can be different on the right and left side of the 
injection electrode, Vd{x) = 6{x)v^—9{x)v^^ and is determined by the corresponding currents 
driven on either side of the injector. For a special case of = v^^ the steady state spin 
density solving Eq. ^ is given by the commonly used expression,^ 
• dt 

Sy{x) = So J ^^^-^ exp[(-x - Vdtf/Wt] eyi^[-t/Ts] x cos{giJ.BBj/h) ; (3) 

Sz{x) is obtained by replacing cosine by sine in the above expression. Assuming the step-like 
discontinuity in the drift velocity at the injection point, which corresponds to our exper- 
imental geometry, the solution of Eq. ^ outside the injection point must have the same 



functional form as the expression (10), up to a normalizing factor. (Outside the injection 
point, Eq. ^ has the same form of a homogeneous differential equation for both the con- 
stant or step-like Vd{x).) The origin of the renormalization due to Vd{x) with a sharp step at 
the injection point is that this form of Vd{x) is equivalent to an additional source/sink term 
in the drift-diffusion equation at the injection point {d9{x)/dx = 5{x)). As confirmed by our 
numerical solution of the drift-diffusion equation, the two normalization factors for the right 
and left spin densities are obtained by matching the spin densities at the injection point 
and by requiring the same total integrated spin density as in the case of the constant drift 
velocity, i.e., dxSy{x) = rsSo/[l + {ubTsY] and dxs^ix) = rsSo{ujBrs) /[I + {ujbTsY]- 
Note that, the conservation of the integrated spin density is valid for spatially independent 
spin-dephasing time and magnetic field in Eq. 



The drift velocities corresponding to our experiments in Figs. 4b, c are given by, 



Id /en A and = {Id + Ib) / enA (see Fig. 4a). Here e is the electron charge, n is the electron 
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density in the GaAs channel, and A is the cross-sectional area of the channel. At the low- 
temperature used in the measurements, the diffusion constant is given by the expression 
for a degenerate semiconductor, D = ji^n/ eg{EF)^ where /Xg is the electron mobility and 
g{EF) is the density of states at the Fermi level in GaAs conduction band with effective 
mass rrf = 0.067. The mobility /ig = 3.5 x 10^ cm^V~^s~^ and density n = 1.1 x 10^'' cm~^, 
and the corresponding diffusion constant D = 2.9 x 10~^ m^s~^ and drift velocities were 
determined using the ordinary Hall measurements in the GaAs channel. (To extract the 
coefficients from the ordinary Hall data we considered A = wt where the width of the 
channel w = 20 /xm and the effective thickness of the conducting GaAs film t = 270 nm.) 

The spin-dephasing time = 1.65 ns is obtained by matching the width of the theoretical 
and experimental Hanle curves. We determined from measurements in the applied out-of- 
plane hard- axis field 5^, i.e., in the geometry where the Overhauser field is negligible. The 
remaining input parameter needed for obtaining the quantitative values of the theoretical 
non-local spin valve Hanle curves, shown in Fig. 4d, is the overall normalization factor of the 
continuous solution of Eq. ^ (or equivalently the value of Sq). This is obtained by matching 
the theoretical and experimental spin densities in GaAs underneath the detection electrode. 
The experimental value is inferred from the difference between the zero field non-local spin 
valve voltages at parallel and antiparallel magnetization configurations of the injection and 
detection Fe electrodes considering ,^1 

AVnl = 27? . (4) 

6e 

Here rj = 0.5 is the spin transmission efficiency of the interface, Ppe = 0.42 is the polarization 
of the Fe electrode, and PcaAs = '^Sy{xd)/n is the polarization in GaAs underneath the Fe 
detection electrode (x = Xd). 

The iSHE is proportional to the z-component of the spin-current given by j'l(x) = 
—DVsz{x) + Vd{x)sz{x). Since j^x) depends on the spatial coordinate we have to con- 
sider also the response function Fcrossi^) of the finite-size Hall cross when interpreting the 
experiments. We performed the numerical evaluation of Fcrossi^) for our sample geometry 
using the conformal mapping theory (see Supplementary information) The measured 
iSHE signal is then proportional to J| = J^^dxj^{x)Fcross{x) / JT^d^Fcrossi^)- The spin 
current and the iSHE voltage are related as, Vh = ewaJ^/a^ where a is the spin Hall an- 
gle and a = nefie is the electrical conductivity of the GaAs channel. The theoretical Vh 
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Figure 4: (a) Schematics of the experimental setup. (b),(c) Experimental non-local spin valve and 
iSHE signals in the in-plane hard- axis field Bx measured at constant spin-injection bias current 
Ib — 300 fiA and at three different drift currents Id depicted in (a). The data were measured 
at the magnetic field sweep rate 3 mT/min. (d),(e) Theoretical calculations of the non-local spin 
valve and iSHE signals. 

plotted in Fig. 4e is obtained by taking a = 1.5 x 10~^ which is a value consistent with 
the estimates of the skew-scattering Hall angle for the disordered weakly spin-orbit coupled 
GaAs channel (see Supplementary information). The value is also consistent with electron 
density dependent spin Hall angles in diffusive GaAs channels reported in optical spin Hall 
measurement s.*^^ Figs. 4b, d and 4c, e demonstrate the agreement we obtain between the 
measured and calculated non-local spin valve and iSHE voltages. The theory successfully 
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describes the dependence of the measured spin signals on both the apphed magnetic field 
and on the applied electrical drift current. 

To conclude, we have demonstrated a spin-transistor device based on the iSHE detec- 
tion of spin currents injected electrically into a semiconductor from a ferromagnetic contact 
and manipulated by an electric field. Our demonstration was made possible by designing 
ultrathin-film Fe electrodes with a strong in-plane magnetocrystalline anisotropy to elimi- 
nate the ordinary Hall signal and by experimentally introducing an electrical spin manipu- 
lation method suitable for diffusive semiconductor channels. We have performed a detailed 
quantitative analysis of the measured iSHE signals based on complementary non-local spin 
valve measurements of the injected spin polarization, and on solving the spin drift-diffusion 
equation and the Hall cross response function. 
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— Supplementary information — 
I. EXPERIMENTAL METHODS 

A. Device fabrication and characterization 

Fig. [5] shows the SEM micrograph of one of the devices used in this work. Electron-beam 
hthography and reactive ion etching were used to pattern the lateral GaAs channel with 
the Hall crosses (grey) and with the magnetic electrodes (white). The injection contacts 
were first defined by Ti/Au using lift-off. Then the Fe/Al double layer was selectively 
etched everywhere else. Next the Hall bar was defined using reactive ion etching. Finally 
the injection contacts were connected with contact pads (outside the figure) using Ti/Au 
air-bridges. 

We performed local tunnelling anisotropic magnetoresistance (TAMR) measurements to 
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Figure 5: SEM micrograph of the device. 



monitor the magnetization orientation of the individual Fe electrodes. Fig. [6](a) shows the 
circuit setup we used for the three-point measurements between an individual Fe electrode 
and the two Au contacts. Current / is sent between the Fe electrode and one Au contact while 
voltage V is measured between the Fe electrode and the opposite Au contact. Resistance 
R = V/I is the resistance of the tunneling contact without the contribution of the channel 
resistance in this setup. 

Fig. [6](b) shows the resistance variation as a function of the external magnetic field 
applied closely to the in-plane hard-axis direction. The black squares in Fig. [6](b) are 
measurement points and the red curve is the theoretical fitting. The simulation is based on 
a single domain model with the total energy Etot/Ms = Be/ 4: * cos^(2(/)) + Bu/2 * cos^ (j) — 
Bext cos((/) — 0^). Ms is the saturation magnetization, Be and Bu are the cubic and uniaxial 
anisotropy fields, receptively, cj) and 0^ are the angles between x (i.e. the[ll0] hard-axis) and 
the in-plane orientation of the magnetization, and the applied magnetic field, respectively. 
Minimizing Etot for a given B determines the position of the magnetization at the particular 
applied magnetic field. Since the magnetization angle dependence of the TAMR follow^ 
sin^((/)), the resistance can be expressed as i? = i?o + Asin^ (j). Rq is the tunneling resistance 
with Fe magnetization along the hard-axis direction, and A is the amplitude of the TAMR 
signal. The good agreement between experimentally obtained values and the simulation 
confirms the single domain behavior of the Fe electrodes. The fitting yields the following 
anisotropy fields: Be = 20 mT and Bjj = 154 mT. These values are consistent with previous 
studies of ultra-thin Fe layers epitaxially grown on GaAs.l^The anisotropy field required to 



16 



n<2> 




-300 -150 150 300 120 130 

B (mT) B (mT) 



Figure 6: (a) TAMR measurement geometry (b) Measured (black squares) and simulated (red 
curve) TAMR signal of an Fe electrode, (c) Switching of two Fe electrodes monitored by the 
TAMR. 

rotate the magnetization from the easy to the hard-axis direction for external field applied 
along the hard-axis is Ba = Bu + 25^ = 194 mT. 

An example of detailed TAMR measurements close to the switching fields for two diflFerent 
Fe electrodes are shown in Fig. [6](c). The switching is seen as the step of the TAMR 
whose position and size depends on the applied field angle and the measured electrode. 
Once the switching fields are known for each electrode, the magnetizations can be switched 
individually by applying a suitable field. 

B. Experimental extraction of the iSHE 

In Fig. [7|(a) we show an example of the extraction of the iSHE signal from the measured 
raw data. The external magnetic field B applied along the in plane magnetic hard-axis is 
swept from positive to negative values at a rate of 3mT/min. The blue curve shown in 
the figure are the measured raw data V^, obtained by setting the magnetization in the Fe 
injection electrode in the positive easy-axis orientation, with a subtracted linear background 
of 0.3/il//mT. The linear component is attributed to the ordinary (Lorentz force) Hall eflFect. 
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The ordinary Hall effect may have two contributions in the geometry of our experiment: One 
due to the vertical part of the trajectory of the electrons passing from the Fe electrode to the 
GaAs channel, and the other one from the lateral part of the transport in the GaAs channel 
due to a small out-of-plane misalignment of the sample and the applied magnetic field. (The 
Hall slope of O.S/iV /mT corresponds to a small misalignment of few degrees, which is within 
the experimental error of orienting the sample for the measurement.) The ordinary Hall 
eflFect may also have a contribution from the stray field of the Fe electrode. Importantly, 
both the contributions from the external field and from the stray field are independent 
of the Fe magnetization orientation (see Fig. [7](b)). Therefore, the ordinary Hall effect 
is eliminated from the measured data by subtracting the measured raw data obtained by 
setting the magnetizations in the Fe injection electrode in the positive and negative easy- axis 
orientations. Consistently the resulting signal, VJP = (V^ — V^)/2, shown as the black curve 
in Fig. [7](a) overlaps with the blue curve. The black and blue curves are of pure spin origin 
and represent the measured iSHE. 

The red curve in Fig. [7](a) shows the antisymmetric part of the iSHE signal, V^^{B) = 
[Vff{B) — Vff{—B)]/2. The antisymmetrization of the signal allows us to correct for the 
residual asymmetry of the absolute value of the signal coming from the fact that measure- 
ments were done at a rate only approaching the equilibrium nuclear spin polarization. When 
coming from the large positive fields, the nuclear polarization is higher compared to its equi- 
librium state, hence it dephases electrons more strongly. When sweeping from zero field to 
the large negative fields, the opposite situation occurs making the absolute value of the 
signal weakly asymmetric with respect to positive and negative fields. By comparing the 
red curve with the blue and black curves we see that the residual asymmetry of the absolute 
value of the signal is relatively weak and is not obscuring our iSHE data. 

Note, that the experiments were performed with a small intentional in-plane misalignment 
(1 — 3°) from the in-plane hard- axis which allowed us to preset the magnetization direction 
without rotating the sample during the measurements. We checked that the obtained iSHE 
signal was independent of this small in-plane misalignment. 
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Figure 7: (a) The blue curve shows the measured raw data Vj, obtained by setting the magne- 
tization in the Fe injection electrode in the positive easy-axis orientation, with a subtracted linear 
background of 0.3/xy/mT. The black curve shows the pure spin signal VJ^ obtained by, = 
(t4 - ^^)/2. The red curve shows the antisymmetrized signal V^%B) = [V^{B) - VJp{-B)]/2. 
(b) Schematics of the stray fields when the magnetization is preset along the positive and negative 
easy-axis orientations and the external hard-axis field is applied. 

C. Determination of the sign of the spin-polarization 

In order to determine the sign of the injected spin polarization with respect to the mag- 
netization orientation of the injection electrode, we analyzed the response of the injected 
electron spins to the nuclear spins which were polarized in a controlled way. This technique 
was previously applied in Ref. 7. The effective magnetic Overhauser field, generated by the 
dynamically polarized nuclear spins (Eq. (1) in the main text) is non-zero when the applied 
magnetic field B is not aligned perpendicular to the electron spin polarization (S), i.e. to 
the magnetization of the injection electrode. At low applied magnetic fields intentionally 
misaligned from the hard-axis direction, the Overhauser field can compensate for the ap- 
plied magnetic field and suppress the resulting spin-decoherence and spin-precession so that 
a satellite peak appears in the Hanle measurement of the non-local spin valve signal when 
the Overhauser and external fields cancel each other. The position of the satellite peak with 
respect to B = determines the sign of the spin-polarization. 

In Fig. [8] we show non-local spin valve Hanle signals measured for the parallel config- 
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Figure 8: (a) Sketch of the nonlocal spin valve Hanle experiment determining the sign of the 
spin polarization: The applied magnetic field B is misaligned with respect to the in-plane magnetic 
hard-axis by ^ 10°. (b) Hanle measurements of the non-local spin valve signal for the parallel 
configuration of injector and detector magnetizations along positive (black) and negative (red) 
easy-axis orientations. B is swept in both cases from +50mT to -50mT at a rate of 3mT/min. 

uration of injector and detector magnetizations along positive (black) and negative (red) 
easy axis orientations. After presetting the magnetizations, a magnetic field misaligned 
from the in-plane magnetic hard-axis by ^ 10° is swept from +50mT to -50mT at a rate 
of 3mT/min. The Overhauser field is positive when the magnetization is oriented along the 
positive easy- axis orientation because the satellite peak appears at S < (black). When the 
magnetization is along the negative orientation, the satellite peak is found at 5 > and the 
corresponding Overhauser field is negative (red). This indicates that the spin polarizations 
of nuclei and injected electrons are opposite to the magnetization orientation of the Fe spin 
injection electrode. 

In the case of the external field applied close to the hard- axis direction, the component 
of S in the direction of the external field is given by the tilt of the magnetization towards 
the applied field. In the case of the opposite sign of accumulated electron spins to the Fe 
magnetization, this results in an enhancement of the effective field acting on electrons as 
introduced in the main text. 
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II. THEORY DISCUSSION 



A. Spin drift-diffusion equation 

We begin with the drift-diffusion equations in one dimension, 

dS,{x,t) ^_I)^Si{x,t) + va{^)Si{x,t)) + +7e,,-,B,-5fc(x,t) = Gi{x,t) (5) 

at Ts 

where D is the diffusion constant, Vd is the drift velocity, is the spin-dephasing time, and 

— * 

7 is the gyromagnetic ratio. G is the injection rate. Using these relations to rescale the 
relevant quantities, x* = T = — , aS* = G* = = ITsB, iT/ = y^'u^, and 

assuming, = Dts the equation can be reduced to: 

^ + V ■ {-VS* + v*,S:) + S: + e,j,B*Sl = ^ (6) 

1. Zero field case 

The zero magnetic field can be solved analytically in a straight forward way. Following 
our experimental set-up let us assume that G{x) — Go5{x)z^ and that there is no magnetic 
field. Then, in the steady-state case (when the time derivative is zero) we have: 

- s;"{x) + K{x)s;{x)y + s;{x) = six) (7) 

Here we have written the equation for the ^/-component of S only, since without magnetic 
field each component is uncouple, and since we are injecting in the ^/-direction, only Sy is of 
interest. Assuming that is constant, 

- S'^ix) + VdS'yix) + Sy{x) = S{x). (8) 

Note that Vd refers to the scaled drift velocity as prescribed before. It is also the only free 
parameter in this differential equation. The solution is S{x) — ^e^^e~^^^\ where a — 
and cj^ = 1 + a^. 

We next consider the case that Vd is given by a discontinuous step function, such that 
'^cz(^) — ^d^{~^) + ^d^{^)' The solution is given by: 

S{x) = — -[e(«"+'^")^^(-a;) + e("''-'^'')^^(x)], (9) 

where = f , = f ,(0;^)' = 1 + (u^f = 1 + (a^. 

21 



2. Non-zero field case 



For the case of finite magnetic field it is more straight forward and transparent to proceed 
in a simple numerical way. For the case of constant drift velocity the solution has the form 

• 1 

Sy{x) = So J ~^J^^ exp[(-x - Vdtf/Wt] exp[-t/r,] x cos{giLiBBj/h), (10) 

with the solution for Sz obtained by replacing cosine by sine in the above expression. A sharp 
step-like form of Vd{x) is equivalent to a source/sink term in the drift-diffusion equation, 
since a discontinuous drift velocity will tend to accumulate spin at the discontinuity. Within 
a constant approximation, the full steady-state solution of the drift-diffusion equation can 
be shown to be normalized to rg5'o/(l + (o^b^sY) the y-component and TsSo{(jJbTs) / {I + 
(ujb^sY) the z-component. Hence, the solution for Sy{x) and Sz{x) assuming a step- 



function behavior of Vd{x) is obtained by Eq. 10 for constant Vd on the left and the right 
of the drift velocity discontinuity with Eq. [TO] multiplied by appropriate constant factors on 
the left and on the right to make Sy{x) and Sz{x) continuous and normalized correctly. 

In Fig. [9] we illustrate the Sy{x) for the situation where the additional drift current {1^ in 
Fig. 4 in the main text) is +100 //A, 0, —100 /xA and the bias current through the injection 
electrode Ib = 300 /xA, as done in the experiment. The spin-current generated by this spin 
accumulation profile is given by 

J, = -DVs,{r) + Vd{f)s^{r). (11) 



B. Hall effect 

The conduction electrons can be modeled through the following effective Hamiltonian: 

H=^ + ydis(r) + AV • (k X VVdis(r)), (12) 

where the m = O.OGTmg and Vdis is the disorder potential modeled by uncorrelated delta 
scatterers of strength Vq and density n^. For GaAs A* = 5.3 K^W^The Hall effect signal 
can be understood within the theory of the anomalous Hall effect ( AHE) . The contributions 
to the AHE in SO coupled systems with non-zero polarization can be classified in two types: 
the first type arises from the SO coupled quasiparticles interacting with the spin-independent 
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Figure 9: Spin polarization profile Py(x) = Sy{x)/n obtained by solving the drift-diffusion equa- 
tions for the experimental parameters of Fig. 4 in the main text. 

disorder and the electric field, and the second type arises from the non-SO coupled part of 
the quasiparticles scattering from the SO coupled disorder potential. The contributions of 
the first type do not dominate the physics of this weakly spin-orbit coupled system, and is 
therefore not included in Eq. [T2j 

The contributions of the second type, i.e. from interactions with the SO coupled part of 
the disorder ,1^^^ are due to the anisotropic scattering, the so called extrinsic skew-scattering, 
and is obtained within the second Born approximation treatment of the collision integral in 
the semiclassical linear transport theoryP^'^ 



a- 



I skew 



xyl 



(13) 



Using the relation for the mobility fi = er/m and the relation between n^, Vq, and r, 
h/r = riiV^m/h^ ^ the extrinsic skew-scattering contribution to the iSHE angle due to a pure 
spin-current, a = Pxy/Pxx ^ (^xy/(^xx-> can be written as 

A*[A>2D[10^'cm-2] 



a 



skew 



2.44 X 10" 



3.8 X 10"^ 



^/i[103cm2/Vs]n2D-i[10iicm-2] 



(14) 



where we have used n2D = nt = 3.0 x IQ-^^cm ^, t = 270 nm is the GaAs film thickness, 
n = 1.1 X 10^'^cm"^, /i = 3.5 X lO^cm^/Vs, and n2D-i ^ 3.0 x lO^^cm"^. 
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C. Hall Response Function 



When dealing with non-uniform currents, it is non-trivial to relate the measured Hall 
voltage signal with the Hall angle or Hall coefficient of the system.l^SEII in geometries where 
the Hall probe width and the channel width are of similar magnitude the current density 
near the the cross can contribute to the Hall signal more significantly and one must solve 
the full equations relating the current density and the fields. For the case where anomalous 
Hall effect is considered in addition to the normal Hall effect we can write 

J(x, y) = p-\-VV{x, y) + J(x, y) x [R^B{x, y) + 47ri?,M(x, y)] , (15) 

where p is the diagonal electrical resisitvity of the layer, B and M are the local magnetic 
induction and magnetization, and i?o and Rg are the normal and anomalous Hall coefficients. 
Here we have assumed a thin film geometry such that = and the problem is reduced to 
two dimensions. 

We then must solve the Maxwell equations for a static magnetic field 

VV = and V- J=0 , (16) 

with the boundary condition that the current is zero at the insulating cross boundaries, i.e. 
j • h = dit the boundaries. 

We follow here a similar procedure as in Ref. 23. Taking /3 = p~^[RoB{x^y) + 
4:7TRsMj_{x^y)] the above equations reduce to 

(1 + f3^))V^V + (1 - /5^) - - 2f3 + - (17) 

\dx dy dy dx J \dx dx dy dy J 

and the boundary condition to V±V = —/3V\\V. Since /3 tends to be small in most systems 
of interest we can treat the problem perturbatively, V = Vq + Vi + . . . ^ whose first two 
components solve 

Wo = 0, (18) 

with dVo/dn = at the boundaries and 

vv; ^ f'^ - f'^, (19) 

ay ox ox oy 

with V±Vi = — /3V//Vo at the boundaries. Solving these equations for the case of a delta- 
like magnetic field at position (x^y) yields the Hall response function which can then be 
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convoluted with the non-constant magnetic field or magnetization to obtain the total Hall 
signal expected. 

In our case, we are considering the response to a pure spin-current which can be considered 
as two fully spin polarized charge currents with opposite polarities and direction. This 
allows us to use the result shown here, only ignoring the small fraction contributing from 
the polarized charge current on the left of the injection point far away, relative to the spin- 
diffusion length, from the detecting Hall bar. 



1. Solution ofVo 

The solution of Vq for the Hall cross bar geometry can be done using the confor- 
mal mapping technique. In here one uses the fact that for any analytical function 
f{z) = u{x^y) + w{x^y) in the complex z-plane, the conjugate functions u and v solve 
the Laplace equation. Then the problem reduces to finding the analytical function that 
solves the boundary conditions of the problem. To do so one can do a conformal mapping 
of the region of interest in the z-plane to a much easier configuration in another complex 
plane, e.g. parallel plate. The conformal mapping preserves the boundary conditions and 
the solution in the z-plane can be obtained by mapping backwards the trivial solution in the 
complex plane. 

For the case of the Hall cross, or any polygon structure for that matter, one uses the 
S wart z- Christ offen transformation. We impose the boundary condition for f[z) — u{x^y) + 
iv{x^y) such that v is equal to zLtt at the boundary edges along the channel. The Swartz- 
Christoffen transformation transforms any interior region of a polygon (even ones with open 
boundaries where the vertex is at infinity) onto the upper half of the complex plane with 
the vortices mapped to points in the real axis. For the present configuration the map reads 



where a, 6, (7, are constants adjusted such that the the vortices map correctly to the right 
places. This maps the problem to a system where the potent ail is ±7r on the real axis which 
can the be mapped simply via a second transformation 
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which is simply the solution of a parallel plate capacitor. Depending on the geometry some 
fraction of the current avoids the central region of the cross bar. 




1.0, b , ^^^^"^^2, C 
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Figure 10: (a) Hall cross response function in the x — y plane for the geometry of the mea- 
sured device, (b) Hall cross response function averaged over the channel width (^/-direction), (c) 
Spin-current profile obtained by solving the drift-diffusion equations for the experimental 

parameters of Fig. 4 in the main text. 

2. Solution of the Hall response function: Vi for a Dirac delta function magnetic field 

Since the equations that we are solving, to first order, retain the principle of superposition 
we can solve for the Hall response function due to a delta function magnetic field. 
The Hall response function is found (see App. in Ref. 23) to be 



The above response function is normalized to 27r. The response function for our experimental 
set-up is shown in Fig. |10| (a). In realistic situations there is no current with about 100 nm of 
the edge and we therefore exclude the sharp part of the response function near these edges. 
Since we are only considering the spin accumulation in one dimension, we can average this 
response function in the y-direction, as shown in Fig. |10|(b). Solving for the spin-current 




(22) 



2 Iw^ _ 0,2 
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from the drift-diffusion equations in one dimension, shown in Fig. uOfc), using Eq. 11 



we then convolute this result with the response function integrated along the ^-direction 
to obtain the measured spin-current J| = J^^dxjl{x)Fcross{x) / JT^d^Fcrossi^) which is 
related to the iSHE voltage by Vh = eujaJl/a as described in the main text. 
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